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Abstract 

In this paper, we consider the perturbed renewal risk process. Systems of integro- 
differential equations for the Gerber-Shiu functions at ruin caused by a claim and 
oscillation are established, respectively. The explicit Laplase transforms of Gerber- 
Shiu functions are obtained, while the closed form expressions for the Gerber-Shiu 
functions are derived when the claim amount distribution is from the rational family. 
Finally, we present numerical examples intended to illustrate the main results. 
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1. Introduction 

Consider a continuous time renewal risk process perturbed by diffusion 

N(t) 

U(t)=u + ct-^2Zi + aB(t), t>0, (1.1) 

i=i 

where u > is the initial capital, c > is the constant rate of premium. The 
ordinary renewal process {N(t),t > 0} denotes the number of claims up to time t, with 

N(t) = max{n >l:V 1 + V 2 -\ h V n < t}. Then V i} i = 1, 2, . . . are the interclaim 

random variables. They are independent and assumed to have common distribution 
function K, density function k, and Laplace transforms k(s) = / °° e~ sx k(x)dx. {Z iy i > 
1} are independent claim-size random variables with common distribution P (such 
that P(0) = 0) and density p. {B(t),t > 0} is a standard Brownian motion with 
B(0) = 0. It is assumed that {N(t)}, {B(t)} and {Zi} are mutually independent and 
that cE(Vj) > E(Zi) providing a positive safety loading factor. 

The perturbed risk model of form (1.1) was firstly introduced by Gerber (1970) 
and has been studied by many authors. See, for example, Dufrense and Gerber (1991), 
Furre and Schmidli (1994), Schmidli (1995), Gerber and Landry (1998), Wang and Wu 
(2000), Tsai and Willmot (2002a, b), Zhang and Wang (2003), Li and Garrido (2005) 
and the references therein. 

Let T = inf{t : U{t) < 0} (with inf{0} = oo) be the time of ruin for risk process 
(1.1). Two important nonnegative random variables in connection with the time of 
ruin T are |C/(T)|, the deficit at the time of ruin, and U(T—), the surplus immediately 
before the time of ruin. Consider a penalty scheme which is defined by a constant w 
if ruin occurs by oscillation and w(U(T—), \ U(T)\) if ruin is caused by a jump. Then, 
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the Gerber-Shiu discounted penalty function at ruin (p{u) is defined by 



<f>(u) = <f> w (u) + w <f> d (u), 



(1.2) 



where for 5 > 0, 



<f> w (u) = E [e~ 5T I(T < oo,U(T) < 0)w(U(T-), \U(T)\)\U(0) = u] , (1.3) 



(with <f> w (0) = 0) is the expected discounted penalty function at ruin caused by a claim, 



(with 0^(0) = 1) is the Laplace transform of the ruin time T due to oscillation. Many 
ruin-related quantities can be analyzed by appropriately choosing special penalty func- 
tion w, for example, let 5 = and w — 1, then 4> w (u) = ip w (u) gives the probability 
of ruin due to a claim and <f>d(u) — ^pd(u) is the probability of ruin that is caused by 
oscillation. 

The evaluation of the Gerber-Shiu discounted penalty function, first introduced in 
Gerber and Shiu(1998), is now one of the main research problemes in ruin theory. See, 
for example, Gerber and Landry (1998), Tsai and Willmot (2002a,b) for the classical 
surplus process perturbed by diffusion, Li and Garrido (2005) for the generalized Erlang 
(n) risk process perturbed by diffusion, Albrecher and Boxma (2005) for the semi- 
Markov model, Willmot (2007) and Landriault and Willmot (2007) for the renewal risk 
model, Lu and Tsai (2007) for the Markov-Modulated process perturbed by diffusion. 

The rest paper is organized as follows. In Section 2 we derive systems of integro- 
differential equations for Gerber-Shiu functions. Section 3 discusses a generalized Lund- 
berg's equation and its roots. And the Gerber-Shiu functions for the model are fully 



and 



<t> d (u) = E[e- 5T I(T < oo,U(T) = 0)\U(0) = u], 



(1.4) 
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analyzed in Section 4. Section 5 contains several numerical examples intended to illus- 
trate the main results. 

2. Integro-differential equations 

More recently, Ren (2007) considered the risk process with phase-type interclaim 
times, i.e., the distribution of the interclaim time K is phase-type with representation 
(a, B, b), where a. and b are row vectors of length n and B is a n x n matrix. That 
is, each of random variables 14, k — 1, 2, . . . corresponds to the time to absorption in a 
terminating continuous-time Markov Chain with n transient states {81,82, • • -8 n } 
and are absorbing state 8q. Let e denote a row vector of length n with all elements 
being one. Then b T = —Be . Following Asmussen (2000), 

K{t) = 1 - ae tB e T , t > 0, 
k(t) = ae tB b T , t > 0, 

and 

poo 

k(s)= e~~ st k(t)dt = a(sl — B) _1 b T . (2.1) 
Jo 

For % — 1, 2, . . . , n, let 0(w; %) denote the Gerber-Shiu function given £7(0) = u and 
Jg 1 - 1 = S h that is, 

(f)(u; i) = E[e- 5T I(T < oo)w(U(T-), \U(T)\)\U(0) = u, J™ =8,], i = 1, 2, . . . , n. 
Then the Gerber-Shiu function may be computed by 

<f>(u) = a 4>{u) 

where 4>{u) = (<f>(u; 1), . . . , <f>(u; n)) T is a column vector of functions. Similarly, we 
write (f> w (u;i) and cf) d (u;i) for the Gerber-Shiu functions at ruin caused by a claim 
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and oscillation respectively, given U(0) = u and = Si. Denoted by (f> w (u) = 
((f) w (u; 1), . . . , (j) w (u; n)Y and d («) = ((j) d (u; 1), . . . , 0d(w; n)) T . 

Our first result gives integro-differential equations for Gerber-Shiu functions. 
Theorem 2.1. Let w > 0. Then, <f> w (u) satisfies the following equation 

r 2 



a -4>'» + c#» + (B — 5 I)0 W («) + 



ck — x)p(x)dx + u(u) 







0(2.2) 



where u;(u) = w(u, x — u)p(x)dx, I = diag(l, 1, . . . , 1), denotes a column vector 
of length n with all elements being 0, and 4> w (0) = 0, while 4> d (u) satisfies 

r 2 



a —<p" d {u) + c4>' d (u) - 5cf> d (u) + B<t> d (u) + 



u 



a (fi d (u — x)p(x)dx 



o 



b T = 0, (2.3) 



with d (O) = e T . 

Proof. Define {J(t),t > 0} by piecing the {J^} together, 

J(t) = {J t {1) }, 0<t<V u J(t) = {.4% i }, V 1 <t<V 1 + V 2 ,... 

Then { J(t)} is Markov. Jacobson (2005) showed that the joint process {(U(t), J(t)),t > 
0} is a Markovian additive process and thus {(U(t), J(t)),t > 0} is a homogeneous 
Markov process. Then we shall use the technique developed in Grandell (1991, P.84) 
and Wu and Wei (2004) to derive the integro-differential equations for (f> w (u;i) and 
4>d{u; i) for i — 1, 2, . . . , n. Let us consider <f> w (u; i) first. Consider a short time interval 
[0, h}. By noting the fact that claim occurs at the time when the state for { J{t), t > 0} 
changes, we separate the three possible cases: 

(1 ) no change of state occurs in [0, h] for the underlying Markov chain {J(t),t > 0} 
and denoted by J[0, h] = Si in this case; 

(2) the state for {J t ,t > 0} changes in [0, h\ but no claim occurs; 

(3) at least one claim occur in [0, h], 



In the case (1). Denoted by the conditional expectation given the initial 

(C/(0), J(0)) = (u,£i) and V(h) = u + ch + aB(h). Let & J and denote the 

natural filtration of processes {</(£)} and {(U(t), </(£))} respectively. For t > 0, let 9 t 
be the shift operators (see, Revuz and Yor (1991, P.34)). It follows from the Markov 
property of both the underlying process {</(£)} and the vector process {(U(t), J(t))} 
that 



E (u,i) [ e -*r /(r < < 0)w(U(T-), \U(T)\)I(J[0,h] = %)] 



El E 



E^[w{U(T-), \U{T)\)e- 5T I{T < oo, U(T) < 0)/(J[0, h] = Si) 



= E< E 



e- Sh I(J[0, h] = £i 



w(U(T-), \U{T)\)e- 5T I{T < oo, U(T) < 0)) o 6 h \& 



h 



I ^ J 



= e~ 5h E{E 



I(J[0,h] = £i) 



E (u(h),J(h)) [ w (u( T _^ \U{T)\)e- ST I(T < oo, U(T) < 0)] \& t 



e- 8h e^ h E[<t> w {V{h)-i)i 



(2.4) 



where bij is the (i, j)th entry of matrix B. 

For case (2) and (3), by the similar argument to that of case (1), we have 



E (u,i) [ e ^ T /(T < oo, f/(T) <0)w(U(T-),\U(T)\)I(J(0,h] ^S u N(h) = 0)] 



= e-^(l-e^) {^-)E[ct> w {V{h)-j)] + o{hl 
^ -On 



(2.5) 
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where o(h)/h — > as h — > 0, and 



[ e -<5T 7 ( T < qo^^) < o) w ([/(T-), |C/(T)|)/(iV(fc) > 1)] 



w(V(/i), x — V(h))p(x)dx 
v(h) . 



+ o(h), (2.6) 



where b{ denotes the ith entry of vector b. 

Summarizing the above analysis, it follows form (2.4), (2.5) and (2.6) that 



<p w {u;i) = {l-5h + b ii h)E[<p w {V{h)]i)} + h bijEl^iVihyj^+hbi^aj 

j=l,j¥=i 3=1 

rV(h) poo 

I 4> w {V{h) -x; j)p(x)dx + / w{V{h), x - V(h))p(x)dx 

JO JV(h) 



E 



+ o(h). (2.7) 



Apply Ito's lemma for jump- diffusion processes (McDonald (2006), Section 20.8.) 
to (f) w {V(h); k), for k — 1, 2, . . . , n we have 

2 

E[<l> w (V(h); k)} = <j> w (u; k) + [c#>; k) + °-$' w {u- k)]h + o(h). 

Substituting the above expressions into (2.7), canceling <f) w (u;i) from both sides, 
dividing h and letting h — > yields a system of integro-differential equations for 4> w (u; i) 
given the initial surplus u and the initial state of the phase- type distribution Sf. 

2 n n 

3=1 3=1 

a j (^J <Pw(u — x; j)p(x)dx + u(u)^j = 0, for i — l,2,...,n. 

Writing the above equations in matrix form we get (2.2) and note that <f> w (0;i) = 
for i = 1,2, ... , n since P(T < oo, U (T) < 0|£/(0) = 0) = 0. Using arguments similar 
to those used in deriving (2.2), it is not difficult to get (2.3) and <pd(0;i) = 1 for 
i = 1, 2, . . . , n since P(T < oo, U(T) = 0\U(0) = 0) = 1. □ 



Remark 2.1. When the distribution of the interclaim time is a generalized Erlang (n) 
distribution, 



/ 



a = (l,0,...,0),B 



-Ai Ai 







\ 



-A 2 A 2 ••• 



V 



For i = 1, 2, . . . , n - 1, from (2.2) that 



-A„ 



/ \ 





.b 1 



• (2-8) 



a 



\i(j) w {u; i + l) = (Xi + $)4> w (u, i) - aj>' w (u; i) - y^(u; i) 



and 



(A n + 5)<p w {u] n) - c(f)' w (u] n) - —^{u] n) = A r 



<t> w {u — x; l)p(x)dx + ui{u) 



which are formulae (7) and (8), respectively, in Li and Garrido (2005). 

Assume that lim u ^ oo e _SM ?i ,(-u) = and lim^oce -5 ™^^) = hold for 9ft(s) > 0. 
Taking Laplace transforms on both sides of equation (2.2) and noting that (f> w (0) = 0, 
we have 



L(^)0Js) = y ^(O)-^( S )b T , see, 



where 



L(s) = ( — s 2 + cs - 5)1 + B + b T ap(s), 

and Lj(s) = e~ su u{u)du. Denoted by Q w (s) = ^-<t>' w (0) — cD(s)b T . Then the vector 
Laplace transforms 4> w (s) can be solved as 



0«(«) = [L( S )]- 1 Q TO ( S ). 



Thus 

= dSlEwI^'^W- seC (2 ' 9) 

where L*(s) is the adjoint of matrix L(s). 

Similarly, assume that lim u ^ oo e~ SM d (-u) = and ]im u _> 00 e~ su <t>' d (u) = hold for 
9fc(s) > 0. We have 

*'W = 5S]EWI L *W*W' seC - (210) 

where Q,(s) = (£#,(<); 1) + + c, £<ft(0; 2) + f a + c, . . . , f #,(0; n) + + c) T . 

It is observed from (2.9) and (2.10) that the explicit expressions for the Gerber-Shiu 
functions are closely related to the roots of equation: det[L(s)] = 0. This is discussed 
in the next section. 

3. Solutions of Lundberg's fundamental equation 

Let a(s) = ^-s 2 + cs — 5. For values of s such that the matrix a(s)I + B is invertible, 
using the same arguments as in Ren (2007), we have 

det[L(s)] = det[a(s)I + B]det[I+ (a(s)I + B) _1 b T ap(s)] 
= det[a(s)I + B][l + a(a(s)I + B)~ 1 b T p(s)] 
= det[a(s)I + B][l-k(-a(s))p(s)}, (3.1) 

where we utilize (2.1) in the last step. Since the matrix a(s)I + B is assumed to be 
invertible, (3.1) indicates that the solutions for det[L(s)] = and the solutions for 
Lundberg's fundamental equation 

a 2 s 2 

k{5-cs —)p(s) = l, sGC, (3.2) 
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as defined in Gerber and Shiu (2005a) and Li and Garrido (2005) are indentical. 
Theorem 3.1. For 5 > 0, Lundberg's fundamental equation in (3.2) has exactly n 
roots, say pi(S, a), p 2 (S,a), . . . , p n (5, a) with a positive real part $t,(pi(5, a)) > 0. 
Proof. The idea of the proof comes from Gerber and Shiu (2005b). Let 7(5) = 
l/k(—a(s)). Then, as pointed out by Gerber and Shiu (2005b), its zero occurs at 
—a(s) = £, where £ ranges over all n eigenvalues of B. However, B is an intensity 
matrix of all transient states in a continuous time Markov chain, it is nonsingular 
matrix with all its eigenvalues having negative real parts (see Corollary 8.2.1 in Rolski 
et al. (1999)). Therefore, we see that 7(5) has exactly n positive zeros and n negative 
zeros. 

Consider a domain that is a half disk centered at 0, lying in the right half of the 
complex plane, and with a sufficiently large radius. For 9?(s) > 0, obviously that 
\p(s)\ < 1. Because 7(s) has exactly n positive zeros. The theorem follows from 
Rouche's theorem if we can show that \j(s) | > 1 on the boundary of such a half disk. 

Since phase-type distribution belongs to the rational family distributions (see Sec- 
tion 4.2., below), 7(5) has the form of the ratios of two polynomials, in which the 
degree of the denominator is less than the degree of nominator. Then (7(5) | > 1 
for \s\ sufficiently large. Now, for s on the imaginary axis, Sft(s) = 0, we have 
|7(s)| > l/\k(-a(s))\ > 1 also. This ends the proof. □ 

Remark 3.1. If 5 — > 0+ then pi(5, a) — > pt(0, a) for 1 < i < n. and Eq. (3.2) becomes 

a 2 s 2 

k(-cs — )p(s) = 1, s e C. 

Evidently is one of the roots. 

In the rest of paper, pi(5, a) are simply denoted by pi for 1 < i < n and 5 > 0. 
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4. Main results 

4-1- Explicit expressions for (j) 

Here, we recall the concept of dividend differences (see Gerber and Shiu (2005a)). 
For a function L(s), its dividend differences, with respects to distinct numbers g±, g 2 , ■ ■ ., 
can be defined recursively as: L(s) = L(qi) + (s — Q\)L[qi, s], L[qi, s] = L[qi, g 2 ] + (s — 

Q2)L[qi, g 2 , s], The definition of the dividend differences obviously can be extended 

to any vector or matrix that is a function of a single variable. For example, for matrix 

Us), 



Hqi,Q2,s] 



Us) - L( gl ) 

S-Qi 

L[gi,s] - L[gi,g2] 

S - Q2 



and so on. 

We assume that the roots pi,p%, ■ ■ ■ , p n are distinct in the sequel. Since <fi w (s) is 
finite for &(s) > 0. From (2.9) that 



Y L*( Pi )<J>' w (0) = L' t (p i )D(p i )b T , for i = 1, 2, . . . , n. 



Therefore 



y L*[ Pl ,p 2 ]^(0) 



L*(pi)o5[pi, p 2 \ + L*[px, p 2 ]uj(p 2 ) 



and recursively 



-L*[ Pl ,p 2 ,...,p n }(f>' w (0) 



5^L*[pi, p 2 , • • • , Pifi[pi, ■ ■ ■ , p r , 



Then we obtain the following results for <p' w (0) and <p' d (0): 
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Theorem 4.1. The differential of Gerber-Shiu functions at zero can be given by 
1 



</4(o) 
<t>'M 



a 2 / 2 
1 

~^/2 



{L*[p!,p 2 ,...,p„]} 1 



^2'L*\p 1 ,p2,...,Pi]w\p i ,...,p n ] 



b\(4.1) 



{L*[P1, P2, ■ ■ ■ , Pn}} 1 { L*[pi, p 2 , • • • , Pn]e T (^-p„ + C) + 



^LV.pa,...,^^}. (4.2) 



Applying the dividend differences repeatedly to the numerator of Eq. (2.9) and 
(2.10), we obtain the following explicit expressions for the Laplace transform of Gerber- 
Shiu functions: 

Theorem 4.2. The Laplace transform of Gerber-Shiu functions are given by 

nr=i( s -Pi) 



4>M = 



det[L(s)] l L * [Pl ' P2 ' " ' ' Pn ' ^ T^ (0) ~ S (^) bT ) " L+ [^' 

rt— 1 

,Pn-l, , P 2,..., Pi \b T ui[ P i,..., Pn ,s] \, (4.3) 

i=i ^ 



and 



ISsjL^ {L*fc, p 2 , ...,„„, .](^(0) + + Oe-) 



<7 



+ y L *[Pi,P2,---,Pn-i,s]e T 



(4.4) 



Proof. By the fact that s = p\ is a root of the numerator in (2.9), we get 

L*(s)Q w (s) = L*( s )Q w ( s )-L*( Pl )Q( Pl ) 

= [L*(s) - L*( Pl )] Q w (s) + L*( Pl ) [Q w (s) - Q„( Pl )] 
= (s - pM L*[pi, s]Q w (s) + L*(pi)Q w [pi, s] 



= (s-pi)^L*^i,s]Q t0 (s) -L*(p!)b T ^[pi,s] L 



(4.5) 
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Further note that s = p 2 is also a root of numerator in (2.9), implying that s = P2 
is a zero of the expression within the braces in (4.5), Then 

L*(s)Q w (s) = (s-p 1 )^[L*\p 1 , S }Q w (s)-L*(p 1 )b T Q\p 1 ,s}] 

- [L*[ Pl , p 2 ]Q w (p 2 ) -L*(pi)b T tD[pi,p 2 ]] 
= {s-p 1 )(s-p 2 )^L*\pi,p 2 ,s]Q w (s)- 

i, . . . , P2, S\ > , 
i=l J 

recursively from the fact s = p 3 , . . . , s = p„_i are roots of the numerator in (2.9) we 
obtain 

n—l f 

V(s)Q w (s) = Y[(s- p i )<L*[p 1 ,...,p n _ 1 ,s]Q w (s) 
i=i ^ 

n— 1 ^ 

i=i J 

Further note that s = p n is a root of the numerator in (2.9), from (4.6) that 

71—1 s 

L*(s)Q w (s) = ]J(s-p i )<L*[p 1 ,...,p n _ 1 ,s}{Q w (s)-Q w {p n )) + 

i=l ^ 

(L*[pi, . . . , p„_i, S] - L*[pi, . . . , Pn-l, Pn])Q«;(Pn) ~ 
n—l v 

X L *[/°l> • • • ' P^i^lPu • • • > Pn-1, S] - 0[pi, • • • , Pn-1, Pn]) f 
i=l > 

n—l s 

= YI(S - Pi) < L*[P1, • • • , Pn-l, s]Q„,[p n , S] + L*[p!, . . . , p n , S] 
i=l ^ 

n—l ^ 

Q™(Pn) -^L" ; [pi,...,p i ]b T cj[p i ,...,p n ,s] L 
1=1 J 

thus formula (4.3) is derived. Formula (4.4) can be proofed in the same way. □ 
4-2. Closed form expressions of <f> for rational family claim-size distribution 

13 



In some cases the functions 4> w (u) and 4> d {u) can be explicitly and analytically 
determined by inversion of (4.3) and (4.4), respectively. Consider the case where the 
claim-size distribution P belongs to the rational family, i.e., its density Laplace trans- 
form is of the form 

r m {s) 

where r m _i(s) is a polynomial of degrees m — 1 or less, while r m (s) is a polynomial of 
degrees m with only negative roots, all have leading coefficient 1 and satisfy r m _i(0) = 
r m (0). This wide class of distributions includes the Erlang, Coxian and phase-type 
distributions, and also the mixtures of these (see Cohen (1982); Tijms (1984)). 
Multiply both numerator and denominator of Eq. (4.3) by r m (s), yielding 

= U M\V( l \ \ r "^ L ^^ ■ ■ ■ ' S K V^(0) - £(Pn)b T ) - r rn {s) 
r m (s)det[L(s)\ { 2 

n— 1 n 

L*[pi, . . . , p n _i, s]b T Q[p n , s] - r m (s) ^ L *[Pi> P2, • • • , Pi]b T £[pi, . . . , p n , s] W4.7) 

i=i > 

Clearly that r m (s)det[L(s)] is a polynomial of degrees m+2n with leading coefficient 

(^-) n . So the equation r m (s)det[L(s)] = has m + 2n roots on the complex plane. By 

Theorem 3.1. and the definition of the rational distribution, r m (s)det[L(s)] = has n 

positive real roots p±, P2, ■ ■ ■ Pn, and m + n negative real roots. Thus we can express 

r m (s)det[L(s)] by all its roots, i.e. 

2 n m+n 

2 



r m ( S )det[L( S )] = C-Y ]J(s - pi) J] (s + R t ), 



i=l i=l 

where all s have positive real parts. For simplicity we assume that these R^s are 
distinct. Cancel the term Yli=i(. s ~ Pi) fr° m both numerator and denominator of (4.7). 
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Consequently Eq. (4.7) can be rewritten as 



?«,(«) = ^ nm i f : (r m (s)L*[^i, pa, . . . , p n , s](y</4(0) - £(p„)b T ) - r m (s) 



n— 1 ^ 

( s ) 5^ L *[^' ^2, • • • , pJb T cD[pi, . . . , p n , s] \ (4.8) 
1=1 J 

It is not difficult to see that the elements in matrix r m (s)L*[pi, p 2 , . . . , p n , s] or 
r m (s)L*[pi, p 2 , . . . , p n _i, s] are polynomials of degrees which are less than m + n. and 
all L*[pi, p 2 , . . . , pi] for i = 1,2, ... ,n are constants. By decomposing the rational 
expressions in Eq. (4.8) into partial fractions, we get 

r m (s)L [p l ,p 2 ,...,p J ,s] _ ^ M\> _ _ 

where for i — 1, 2, . . . , m + n are coefficient matrices with 

M M = ^(-W..^- forj = „_ 1>n> (4 . 9) 



and 



/ \ ra+n „ 



where Gj for i = 1, 2, . . . , m + n are coefficients given by 



r m (-Rj) 



Thus, by partial fraction Eq. (4.8) can be expressed as 

-, m+n |- 

0„(«) = pT^T E JT^{ M ^^^) " M J n_1) b T -[Pn, s] 

n—l \ 

- Gi^Vlp!, p2, . . . , pt]b T uj\p e , ...,p n ,s] > (4.11) 
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and similarly, 



In order to determine the explicit Laplace inverse of 4> w {s) and (f>d(s), we define 
the same operator T r as in Dickson and Hipp (2001), i.e., for an integrable real valued 
function / with respect to a complex number r (§ft(r) > 0): 

POO 

TJ(x) = / e- T <*-*)f(y)dy, x>0. 



It is clear that the Laplace transform of /, f(s), can be expressed as T s /(0), and for 
distinct ri, r 2 G C 

T T1 T T J(x) = T T2 T T1 f(x) = T * J{X) " rrJ(x) , a;>0. 

r 2 — ri 

Properties of this operator can be found in Li and Garrido (2004). Gerber and 
Shiu (2005a) also presented the following useful result on the relationship between the 
operator T T and the corresponding dividend difference: 



(t[Tr t )f 



i=l 



(0) = (-ir +1 /[r 1 ,r 2 ,...,r n ]. 



(4.13) 



By (4.13) and T s f(0) = f(s), Reserving the Laplace transform of Eq. (4.11) and 
Eq. (4.12), we get the following theorem. 

Theorem 4.3. If the claim-size distribution belongs to the rational family, the Gerber- 
Shiu functions are given by: 



m+n 



n— 1 



(-ir- e L*\ Pl ,p 2 ,...,p e ]b T (l[T pj )u;(u) 



(4.14) 
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where * is the convolution operator, the constants (j = n — 1, n), Gi and Q w (p n ) = 
^4>' w (0) - Q(p n )b T are obtained by Eq. (4.9), (4.10) and (4.1) respectively. And 

i m+n , 2 } 

4>M = E e-^\M^Q d ( Pn ) + ^Mt% T j, (4.15) 

where Qd(p n ) = ^r0d(O) + (^rPn + c)e T can be calculated from (4.2). 

4-3. Ruin probability 

This subsection illustrates the application of the previous results in a special case 
that 5 = 0, w = 1, w(x,y) = 1, p(x) = fie~^ x for x > and the interclaim times 
follow (2.8) with n=2. Then (4.14) and (4.15) become the ruin probabilities caused by 
a claim and oscillation respectively. We now have pis) = (3/(s + (3), oo(u) = e^ u for 
u > 0. The matrix 



Us) 



( a2 2 \ 
^-s -|- cs — Xi X\ 

@>v_ a 2 2 I _„ _ \ 



has exactly two positive real roots p\ = 0, p2 and three negative real roots —Ri, —R2 
and -R 3 . It follows form (4.1) and (4.2) that 

2 

QM = y</4(o)-£(p 2 )b T 

P2 + C — P2 + C-— , (4.16) 



/3(/3 + p 2 )(^p 2 + c) 2 V2^ ' * 2^'^ /3 + p 2 



2 

Q d (P2) = y&(0) + (y P 2 + C)b 

Ai + A 2 



p 2 + 2c/a2\ v 1 1 ( / 3 + p2 )(^ p2 + c )J • (417) 
Moreover, noting that ip w {u) = otif) w (u) = if) w (u; 1) and ipdiu) = ipd(u; 1), together 
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with (4.14) and (4.15) give the following formulae for ruin probabilities 

3 



1 



E 



A1A2 



o*/2)* ^ (3((3 + p 2 )U i¥i (Ri -Ri)V fP2 + c 



1 , T(P~ R i)\--R* 



y-R t u _ e -f3u 



P-Ri 



1 Yl^i( R e - Ri) I ^p2 + c 



Xl + X2 -R l + - 9 )e- R ^ 



(4.18) 
(4.19) 



5. Numerical Examples 

In this section, we will present some numerical examples. In all calculations c = 
1, a — 1, wo — l,p(x) = e~ x for x > are fixed. Let the interclaim times be distributed 
with phase-type representation 

/ -1 ±\ 



a = (1,0), B 



V 







/ 



V 4 / 



Then the mean E[V\] = 9/8, which indicates the relative safety loading is |. 
Example 5.1. (Ruin probability) Let 5 = 0,w(x,y) = 1 for x > 0, y > 0. In this 

case, we have 



L(s) = 



\ 



1 s-1- 

4 

s+l 



1 



2(s+l) 



\ 



4f + s-4 



the Lundberg's equation det[L(s)]=0 has roots: 



pi = 0, p 2 = 2.06412, 

and i?x = 3.90909, R 2 = 3.0744, i? 3 = 0.0806231. 
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Then from (4.14) and (4.15) the different ruin probability components become: 



4> w (u) = 0.27603e- 3 - 90909u - 0.8912 e - 3 - 0744 " + 0.6151 e -' 




0.1 



i 1 1 1 1 1 1 1 1 

2 4 6 8 10 12 14 16 

u 



Fig. 1 . Decomposition of the ruin probability. 



With the help of Matlab, we get Fig. I. for these ruin probabilities for different 
values of u, as well as their decomposition into the ruin probabilities due to claims and 
those due to oscillations. From the graph it can be observed that the ruin probability 
due to oscillations is a strictly decreasing function (from 1 to 0) of the initial surplus 
u. Moreover, when u is small, it decreases sharply, while it decreases slowly when u 
is large. By contrast, the ruin probability due to claims increases quickly at first but 

19 



then deceases slowly after that. 

Example 5.2. (The Laplace transform of the ruin time) When w — 1, we give Fig. 2. 
of E\e- 5T I(T < oo)|C/(0) = u] in the case 5 = 0.1 and 5 = 0.2. 

1.4 
1.2 
1 

(— 

| 0.8 

£= 
CD 

8 0.6 

CO 

0.4 
0.2 


2 4 6 8 10 12 

u 

Fig. 2. The Laplace transform of T. 

As expected, the Laplace transform of the ruin time is high for low 5, the force of 
interest. 
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